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Abstract 

A model-independent search for three-jet hadronic resonance production in pp col- 
lisions at a center-of-mass energy of 7 TeV has been conducted by the CMS Collabo- 
ration at the LHC, using a data sample corresponding to an integrated luminosity 
of 35 pb^^. Events with high jet multiplicity and a large scalar sum of jet trans- 
verse momenta are analyzed. The number of expected standard model background 
events is found to be in good agreement with the observed events. Limits are set on a 
model describing the production of R-parity-violating supersymmetric gluino pairs, 
and gluino masses in the range of 200 to 280 GeV/c^ are excluded at a 95% confidence 
level for the first time. 
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Studies of multijet final states, although experimentally challenging at hadron colliders, are 
sensitive to many extensions of the standard model (SM). For example, variations of techni- 
color models, resulting in heavy colored fermions that transform as octets under SLr(3)c, have 
been proposed in a variety of forms [1 -4 J. Other models incorporate R-parity-violating (RPV) 
decays of supersymmetric gluinos to three-quark final states, where the gluino represents a 
colored adjoint Majorana fermion ||5H7|. In all cases, these high mass resonances can be pair- 
produced, yielding a six-jet final state. Recent results from the Tevatron provide limits on 
gluino RPV decays for masses below 144 GeV/c^ (SI. 

This Letter presents the first results of a dedicated search for three-jet hadronic resonances 
in multijet events in pp collisions. The results are based on a data sample of proton-proton 
collisions at ^/s = 7 TeV, corresponding to an integrated luminosity of 35.1 ± 1.4 pb~^ 191, 
collected with the Compact Muon Solenoid (CMS) detector [10] at the Large Hadron Collider 
(LHC) in the running period from March through November 2010. Events with at least six jets, 
each w^ith high transverse momentum (pt), are selected and investigated for evidence of three- 
jet resonances consistent with strongly coupled supersymmetric particle decays. The event 
selection criteria are optimized in the context of the gluino model mentioned above. However, 
the generic features of the selection criteria provide a robust model-independent basis that can 
be applied to many extensions of the SM. 

The CMS detector is a multi-purpose apparatus, described in detail in Ref. [10 1. Here, we briefly 
describe the subdetectors most relevant to this analysis. The high-resolution silicon pixel and 
strip tracker provides charged tracking coverage for |/y| < 2.4, where t] = — ln[tan(0/2)] is 
the pseudorapidity and 9 is the polar angle measured with respect to the counterclockwise 
proton beam direction. Immersed in the 3.8 T magnetic field of the superconducting solenoid, 
the tracker provides transverse momentum resolution of approximately 1.5% for charged par- 
ticles with Pt ~ 100 GeV/c. Energy deposits of the jets are measured using electromagnetic 
(ECAL) and hadronic (HCAL) calorimeters. The ECAL has a barrel part and two endcaps, 
is composed of finely segmented crystals and has an energy resolution of better than 0.5 % 
for unconverted photons w^ith transverse energies above 100 GeV. The ECAL barrel covers the 
pseudorapidity range \rj\ < 1.4 with a granularity of Arj x Acp = 0.0174 x 0.0174, w^here (p 
is the azimuthal angle, and the endcaps cover 1.4 < \rj\ < 3.0 w^ith a granularity that de- 
creases to 0.05 X 0.05 for \r]\ ^ 3.0. A preshower detector consisting of two planes of sili- 
con sensors interleaved with a total of three radiation lengths of lead is located in front of 
the ECAL endcaps. The HCAL extends up to |?7| ~ 5.0 and its central and endcap regions 
consist of brass /scintillator sampling calorimeters that cover |//| < 3.0 with a granularity of 
Ari X A^ = 0.087 x 0.087 for central rapidities. The energy of charged pions and other quasi- 
stable hadrons can be measured with the calorimeters (ECAL and HCAL combined) with a 
resolution of A E/£ ^ 100 % y/E [GeV] © 5%. 

Events are recorded using a two-tier trigger system. Objects satisfying the requirements at 
the first level (LI) are passed to the High Level Trigger (HLT) where the total recorded rate is 
limited to about ~350 Hz. Triggers based on the sum of all transverse energy from jets (Hj), 
reconstructed with only calorimeter information, are used to select recorded events. For the LI 
trigger, the Hj threshold is 50 GeV. The corresponding threshold for the HLT varies between 
100 and 150 GeV, depending on the run period. 

The CMS particle-flow algorithm fill uses calorimeter information and combines it with re- 
constructed tracks to identify individual particles such as photons, leptons, and both neutral 
and charged hadrons within the jets. The energy of photons is directly obtained from the cal- 
ibrated ECAL measurement. The energy of electrons is determined from a combination of the 



track momentum at the main interaction vertex, the corresponding ECAL cluster energy, and 
the energy sum of all bremsstrahlung photons attached to the track. The energy of muons is 
obtained from the corresponding track momentum. The energy of charged hadrons is deter- 
mined from a combination of the track momentum and the corresponding ECAL and HCAL 
energy, corrected for zero-suppression effects, and calibrated for the nonlinear response of the 
calorimeters. Finally, the energy of neutral hadrons is obtained from the corresponding cali- 
brated ECAL and HCAL energy. The particle-flow objects serve as input for jet reconstruction, 
performed using the anti-A;^ algorithm IIT2II with a distance parameter of 0.5 in rj-cp space. 

Jet energy scale corrections IIT3l derived from Monte Carlo (MC) simulation are applied to ac- 
count for the nonlinear and nonuniform response of the calorimeters. In data, a small residual 
correction factor is included to correct for differences in jet response between data and simu- 
lation. The combined corrections are on the order of 5-10%, and their corresponding uncer- 
tainties range from 3 to 5%, depending on the measured jet's pseudorapidity and energy. Jet 
quality criteria llMll are applied to remove misidentified jets arising primarily from calorimeter 
noise. For both data and simulated signal events, more than 99.8% of all selected jets satisfy 
these criteria. 

Pair-produced gluinos are used to model the signal. Gluino production and decay are sim- 
ulated using the PYTHIA I.15J MC program (v6.420), where each gluino decays to three jets 
through the A^^s quark RPV coupling. This coupling is set such that the branching ratio B 
of the gluino to three light jets is 100%. The mass of the gluino is varied between 200 and 
500 GeV/c^ in 50 GeV/c^ steps. The leading-order cross section from pythia is 325 pb for a 
gluino mass of 200 GeV/c^, falling to ~1 pb for a gluino mass of 500 GeV/c^. For the genera- 
tion of this signal all superpartners except the gluino are taken to be decoupled [7J, the natural 
width of the gluino resonance is taken to be much smaller than the resolution of the detec- 
tor, and no intermediate particles are produced in the gluino decay. The next-to-leading-order 
(NLO) correction factors (K factors), w^ith values ranging from 1.7 to 2.2, are calculated using 
the PROSPINO [161 program and are applied to the leading-order cross sections. Simulation of 
the CMS detector is performed using Geant4 [17J. 

Events recorded w^ith the Hj trigger are required offline to have a good reconstructed primary 
event vertex |18|. Pair-produced three-jet resonances naturally yield events with high jet mul- 
tiplicity and large transverse energy. Thus we require events to contain at least six jets, and that 
the total scalar sum of the pr of those jets is larger than 425 GeV/c. The latter requirement en- 
sures that the trigger is fully efficient for these events. Jets are required to have pr > 45 GeV/c 
and \t]\ < 3.0, which also minimizes the effects from multiple proton-proton interactions. 

To reconstruct the gluinos, the six highest-pr jets are combined into all possible unique triplet 
combinations, resulting in 20 combinations of jet triplets. For signal events, each of the pair- 
produced gluinos corresponds to one of these 20 jet triplets, even in the case w^here all six 
jets come solely from the decay of these particles, leaving the 18 uncorrected jet triplets as 
combinatorial background. Thus, the overall background arises not only from SM events, de- 
scribed by quantum chromodynamics (QCD), but also from spurious jet triplet combinations 
in signal events themselves. We impose additional requirements on each triplet to increase the 
signal sensitivity, w^hile retaining as many signal triplets as possible. The invariant mass of 
background triplets is found to scale with the respective scalar sum of jet pr, while for signal 
triplets the mass is constant. To reduce background, we therefore require each jet triplet to 
satisfy the follow^ing relation: 

3 
%;■ < ElPrI,- - A, (1) 

i=l 



1200 



CMS Simulation \fs = 7TeV 




■~\-W-^'- / 



/ 



/ 



y^ 



y 



-.y 

-^250 GeV/c^ Gluino Model 
20 Triplet Combinations 



.S>10 

c 
m 




10^ 



10^ 



200 



400 



600 



800 



1000 1200 1400 

Jet Triplet Scalar p^ (GeV/c) 



Figure 1: Simulated triplet jet invariant mass Myy^ versus the triplet scalar pj- of all 20 triplets, 
for a gluino mass of 250 GeV/c^. All triplets falling to the right of the red dashed line pass the 
requirement of Eq. IT] In the insert, the invariant mass distribution for the same gluino mass is 
shown both before and after Eq. ITjis imposed. 

where M^yy is the triplet invariant mass, Y^i^\ \VT\i is the scalar sum of jet pr in the triplet (triplet 
scalar pj-), and A is an offset adjusted to optimize signal sensitivity. Figure [T] shows the sim- 
ulated triplet invariant mass versus the triplet scalar p^ for a gluino mass of 250 GeV/c^, and 
the insert displays the invariant mass distribution before and after the requirement. For each 
event, all 20 triplet combinations are included. The value of A is determined by maximizing 
the ratio of the number of signal triplets to the sum of the number of signal plus background 
triplets in a one standard deviation {&) window around the center of the gluino mass peak. A 
common value of A = 130 GeV/c^ is taken for all gluino masses considered, which gives an 
efficiency in signal events for triplets of 1 to 5%, and a background triplet selection efficiency 
of less than 0.05%. 

Even after the final selection, background remains from both QCD multijet events and uncor- 
related triplets in gluino signal events. The latter only contribute minimally, and the shape of 
their distribution is found to be consistent with that of the dominant background, from QCD 
multijet events. These QCD multijet events arise from hard two-particle interactions combined 
with initial- and final-state radiation in the form of gluon jets. Although the cross section falls 
with increasing jet multiplicity (Njet), the underlying kinematic distributions are essentially the 
same among these events. Thus, we use a reseated mass distribution of triplets in events with 
Njet = 4, where the signal contributions are minimal, to estimate the shape of the background. 
Specifically, we select events in data with Njet = 4 that satisfy all other selection criteria, form 
jet triplets, and require each to pass Eq. ITj The M^yy values of these triplets are multiplied by 
the ratio of the average triplet scalar pj in data for events with Njet > 6 to the events with 
Njet = 4, to account for expected minor kinematic differences between the two samples. The 
resulting Myyy distribution is then fit to an exponential function of the form: e^'^^^^^'w, where 
Pq and Pi are free parameters. The slope Pi of the exponential function in the Njet ^ 6 sample 



is constrained to be equal to that found for the scaled Njet = 4 fit within its uncertainties. This 
is verified in QCD simulation, and as a cross-check in data, we apply this procedure to predict 
the shape of the Mjjj distribution for an Njet = 5 sample, where the QCD multijet background 
is also expected to dominate, and find good agreement. To verify that the choice of the back- 
ground model does not bias the derived limit, the exponential function is tested on an Njet > 6 
sample, defined by the standard selection criteria without the requirement of Eq. IT] imposed. 
The parameterization is found to be in agreement with the data in the fitted region, with the 
slope of the fit consistent with those of the Njet = 4 and Njet = 5 samples. 

To estimate the number of signal events expected after all selection criteria are applied, the 
sum of a Gaussian function that represents the signal and the exponential function that mod- 
els the background is fitted to the simulated Mjjj distribution for each gluino mass. The fit is 
performed in the range 170 < Mjjj < 800 GeV/c^. The width of the Gaussian function mod- 
eling the signal varies according to the detector resolution, and gluinos of mass from 200 to 
500 GeV/c^ correspond to widths from 10 to 25 GeV/c^. The integral of the Gaussian compo- 
nent provides the estimate for the expected number of signal triplets produced, and the value 
of this integral, divided by the number of signal events generated, determines the signal ac- 
ceptance for each gluino mass. The signal acceptance is parameterized using a second-order 
polynomial as a function of gluino mass, and the acceptance ranges from 0.4 to 5% as the gluino 
mass increases from 200 to 500 GeV/c^. 

The systematic uncertainty on the signal acceptance is evaluated in the following way. An un- 
certainty related to the jet energy scale [13] is evaluated by varying the jet energy scale correc- 
tion within its uncertainties, then recalculating the acceptance for different gluino mass values. 
The largest difference with respect to the nominal acceptance is taken as the systematic uncer- 
tainty and ranges from 7 to 16%. The amount of initial- and final-state radiation is increased 
and decreased according to a standard prescription |19|, and the associated uncertainty is quan- 
tified in a similar manner to that described for the jet energy scale uncertainties. The difference 
of 2 to 4% with respect to the nominal acceptance is taken as the systematic uncertainty. To 
determine the effects of additional proton-proton interactions on the signal acceptance, signal 
samples are generated with the number of interactions per bunch crossing in the simulation 
set to the average of their distribution in the data. Applying the acceptance calculation on this 
sample leads to differences of 1 to 6%, which are taken as uncertainties. These contributions, 
combined with those from the luminosity measurement (4%) and choice of parton distribu- 
tion function set (4%), give a total systematic uncertainty on the signal yield between 10 and 
19%, depending on the value of the gluino mass. Other effects, such as additional background 
parameterizations and variations of the fit range, were also tested and found to be negligible. 

Figure [2] shows the three-jet invariant mass distribution for the Njet ^ 6 sample w^ith all selec- 
tion criteria applied, and the exponential fit superimposed. The simulated signal distribution 
for a gluino mass of 250 GeV/c^, normalized to the integrated luminosity of the data sample, is 
also shown. Because agreement is observed between the data and expected QCD background, 
a limit-setting procedure is performed. 

Upper limits are placed on the cross section as for the production of three-jet resonances in 
the Njet > 6 sample using a Bayesian approach. The background model parameters and their 
corresponding uncertainties are taken from the fit of the exponential function to the Njet > 6 
distribution, constrained by the Njet = 4 sample, with all selection criteria applied. The un- 
certainties on the two parameters that describe the background shape, namely the exponential 
slope and normalization, are included as Gaussian priors. The central value is set to the best fit 
value and the w^idth to one standard deviation. The range is truncated at ±3(7". In addition to 
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Figure 2: Three-jet invariant mass distribution of triplets passing all selection criteria for the 
Njet > 6 data sample. An exponential function representing the background shape, constrained 
from the Njet = 4 distribution, and the expectation for the 250 GeV/c^ gluino signal are also 
shown. 



Table 1: Observed and expected 95% CL upper limits on the cross section times branching ratio 
for the pair production of gluinos with masses (Mjjj) ranging from 200 to 500 GeV/c^. 
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the background parameters, priors are included for the acceptance and integrated luminosity. 

The integrated luminosity, acceptance, signal width and the two parameters of the exponential 
background distribution are all treated as nuisance parameters. The likelihood is combined 
with the prior and nuisance parameters, and then marginalized to give the posterior density 
for as- Integrating the posterior density to 0.95 of the total gives the 95% confidence level (CL) 
limit for ^75. Marginalization and integration of the posterior density are performed w^ith a 
Markov chain MC integration technique using RooStats 12011 . 

To determine the expected limits, a large set of pseudo-experiments (PEs) is generated using 
the background-only model. For every PE, each of the two parameters associated with the 
exponential is varied by generating a random number distributed according to a Gaussian 
probability distribution function centered at the central value, with a width corresponding to 
the associated uncertainty. The total number of events in a given PE is extracted according 
to the Poisson distribution w^ith mean value equal to the number of events predicted by the 
exponential function in the fitted range. The same upper limit calculation performed on data 
is repeated for each PE at each mass, and the median of the upper limit distribution for all PEs 
is the expected limit. 

The observed and expected 95% CL upper limits on the gluino pair production cross section 
times branching ratio as a function of gluino mass are presented in Fig. Island Table [ij The cor- 
responding 95% CL lower limit on the gluino mass is set by finding the mass value at which the 
95% CL limit line crosses that of the NLO gluino cross section. We thus exclude at the 95% CL 
gluino masses in the range 200 to 280 GeV/c^, with an expected lower limit of 270 GeV/c^. The 
most significant excess occurs for a mass around 390 GeV/c^, corresponding to a significance 
of 1.9 standard deviations, when the so-called look-elsew^here effect ||2T) is taken into account. 

In summary, a search for three-jet hadronic resonance production in pp collisions at a center- 
of-mass energy of 7 TeV has been conducted by the CMS Collaboration, using a data sample 
corresponding to 35 pb^^. Events having the properties of high jet multiplicity and large scalar 
sum of jet Pt, which are expected signatures of high mass hadronic resonances, are analyzed for 
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Figure 3: Observed and expected 95% CL upper limits on the cross section for gluino pair 
production through RPV decays, where the branching ratio of the gluino to three jets is 100%. 
Also shown are the ±lcr and ±2cr bands on the expected limit, as well as the theoretical NLO 
cross section for gluino production. 



the presence of signal events with a model-independent approach. The number of expected SM 
background events is found to be in good agreement with the observed events. The production 
of gluinos decaying through the A^^s RPV coupling is excluded for masses between 200 and 
280 GeV/c^ at 95% CL. These are complementary to recent results from the Tevatron, which 
rule out gluino masses below 144 GeV/c^ ||8l, and these limits are the first from a dedicated 
search of this kind in pp collisions. 
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